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ABSTRACT 

Aims. We aim to describe the pre-main sequence and main-sequence evolution of X-ray and extreme-ultaviolet radiation of a solar 
mass star based on its rotational evolution starting with a realistic range of initial rotation rates. 

Methods. We derive evolutionary tracks of X-ray radiation based on a rotational evolution model for solar mass stars and the rotation- 
activity relation. We compare these tracks to X-ray luminosity distributions of stars in clusters with different ages. 

Results. We find agreement between the evolutionary tracks derived from rotation and the X-ray luminosity distributions from obser¬ 
vations. Depending on the initial rotation rate, a star might remain at the X-ray saturation level for very different time periods, from 
» 10 Myr to »:300 Myr for slow and fast rotators, respectively. 

Conclusions. Rotational evolution with a spread of initial conditions leads to a particularly wide distribution of possible X-ray 
luminosities in the age range of 20-500 Myrs, before rotational convergence and therefore X-ray luminosity convergence sets in. This 
age range is crucial for the evolution of young planetary atmospheres and may thus lead to very different planetary evolution histories. 
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1. Introduction 


High-energy radiation from solar-like main-sequence (MS) stars 
decays in time owing to stellar spin-down. The early Sun’s X- 
ray (a; 1 - 100 A) and extreme-ultraviolet (« 100 - 900 A) 
emissions could thus hav e exceeded the pre sent-day Sun’s level 
by orders of magnitude dRibas et al.ll200-5ll . By driving atmo¬ 
spheric erosion, such extreme radiation levels were critically 
important for both t he primordial hydrogen atmospheres (e.g., 
LammereialJ 2014 ) and the secondary nitrogen atmospheres 
iLichtenegger et all 2010l) of solar system planets. As a conse¬ 
quence of higher solar activity levels, stronger winds would have 
added to atmospheric mass loss through non- t hermal processes 
such as ion pick-up (iKislvakova et al.l 120131 : iKislvakova et al.l 
2 ^ . 


Magnetic activity is strongly coupled to rotation via a stellar 
dynamo, such that the total X-ray luminosity decays with in¬ 
creasing rotation period, P, as Lx oc P^^ to oc P~^; for small P, 
Lx saturates at Lx ~ lO^^Lhni (Lhn i being the stellar bolometric 
luminosity; see lWright et aljl2011 =Wll). Since old er (> 1 Gyr) 
stars spin down in time approxi mately as P oc dSkumanichl 
Il972h . Lx decays as Lx t ' ^ dGiidel et al.l 1 19971) . This evo¬ 
lutionary trend has commonly been formulated using regression 
fits to Lx of stars with known ages, typically starting at the sat- 
urati on level close to the zero-age main- sequence (ZAMS; see, 
e.g.. lGudel et al.ll 19971:[Ribas et al.ll200-5h . 


However, stars in young stellar clusters have a wide 
spread in rotation rates, Q, in particular at ages younger than 
500 Myr before they co nverge to a unique mass-dependent value 
dSoderblom et al.iri993h . As a consequence. Lx values also scat- 


ter ov er a wide range among such stars (e.g., IStauffer et al.l 
Il994h . and the age at which stars fall out of saturation de¬ 
pends on the star’s initial Q. Given the importance of high- 
energy radiation in this age range for planetary atmosphere evo¬ 
lution, a unique description with a single radiation decay law 
is problematic and needs to be replaced by a description of the 
L\ distribution and i ts long-term evolution dPenz et al.l 120081 : 
iJohnstone et al.llT015ll . spanning a wide range of possible evo¬ 
lutionary tracks for stars with different initial fi. 

In this Letter, we use a rotational evolution model to predict 
such luminosity distributions as a function of age, based on a 
range of initial Q, and we show that these predictions agree with 
the observed time-dependent scatter of Lx- We derive a radia¬ 
tive evolution model based on the full range of rotation histories 
for a solar-mass star, and thus hnd a description of the possi¬ 
ble past histories of our own Sun, which is useful to model the 
corresponding evolution of so lar-system planetary at mospheres. 
This Letter is an extension of IJohnstone et al.l (1201 5l) . who sim¬ 
ilarly estimated eyolutionary tracks for wind properties. In this 
Letter, we concentrate mostly on 1 Mq stars, and will extend this 
to other stellar masses in future work. 


2. Rotation and Radiation Models 

As in preyious studies (e.g. iGallet & Bouyieill2Ql^ . we con¬ 
strain our rotation models by assuming that the percentiles of 
the rotational distributions for star clusters with different ages 
can be combined to estimate the time eyolution of a star’s ro¬ 
tation r ate. We consider only s tars in the mass range 0.9 Mq to 
1.1 M,7^. l.lohnstone et ^ (l20L5h collected measured rotation pe- 
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Fig. 1. Comparisons between observed and predicted distributions of Lx at ages of 150 Myr { left) and 620 Myr (righ t). The grey symbols show 
predicted distributions calculated using Eqn.[T]and the distributions of rotation rates derived bv l.Tohnstone et al.1 (l201.5h . The red circles and green 
triangles show detections and upper limits for stars in the Pleiades (left) and Hyades (right). The horizontal blue lines show detection thresholds. 
The upper line of stars in the theoretical distributions is caused by the stars whose rotation rates lie above the mass dependent saturation threshold. 


riods of over 2000 stars in clusters of ages 150, 550, and 1000 
Myr on the MS, giving observational constraints on the per¬ 
centiles at these ages (with 230,134, and 36 stars, respectively at 
the considered ages). We use additional constraints for pre-main 
sequence (PMS) rotation from the «i2 Myr cluster NGC 6530 
(28 stars: iHenderson & Stassun 201^ and the » 12 Myr cluster 
h Per (117 stars: iMoraux et al.ll20i3r For NGC 6530, the 10th, 
50th, and 90th percentiles are at 2.1Q.Q, 6.2fio, and 35.1Go, re¬ 
spectively (assuming Qq = 2.9 x 10“® rad s“'), and for h Per, 
they are at 3.4^0, 8.4f20, and 76.OG©, respectively. 


We use an extens ion of the rotational evolution model of 
iJohnstone et al. (12015 ). For th e wind torque, we use the formula 
derived bv iMatt et al. (1201 2h which relates the wind torque to 
stellar parameters, the star’s dipole held strength, Baip, and the 


urate at a Rossby number of Ro = 0.13, as suggested by the sat¬ 
uration of X-ray emission (Wll), where Ro - Pmtlf* and r* 
is the convective turnover time. For M*, we use the scaling law 
derived by Johnstone et al. (2015), which is derived by htting the 
rotational evolution model to observational constraints. We mod¬ 
ify the scaling law by relating M* to Ro; this allows us to take 
into account the change in t* on the PMS. Since we only con¬ 
sider solar mass stars, the scaling relation derived by Johnstone 
et al. (2015) can be rewritten as M* cc R^Ro ‘‘. We hnd a - 2 
provides a good ht to the observational constraints (which is 
larger than the value of a = 1.33 found by J ohnstone et al. 2015) . 
For Bdip, we use the scaling law derived bv IVidotto et al.l (l2014h 
of Bdip oc Ro 


- 1.32 


To reproduce the spin up on the PMS due to the de¬ 
crease in the stellar moment of inertia, previous studies 
have found that core-enve lope decoupling must be included 
(iKrishnamurthiet n]ll997h. We use the cor e -envel ope decou¬ 
pling model described by iGallet & Bouvierl (l2015h and adopt 
coupling timescales of 30 Myr, 20 Myr, and 10 Myr for the 10th, 
50th and 90th percentile tracks, respectively, which we hnd give 
us good agreement between the rotational evolution model and 
the observations. Finally, we assume that during the hrst few mil¬ 
lion years, the stellar rotation rates do not evolve with time due to 
‘disk-locking’, i.e. magnetic interactions with the circumstellar 


disk. We assume disk-locking timescales of 10 Myr, 5 Myr, and 
2 Myr for the 10th, 50th and 90th percentile tracks, respectively. 

To predict Lx along the rotation tracks, we use the relation 
derived from MS stars by W11, 


Rx = 


CRo^, if Ro > Rosat, 
J^X.sat. if Pt) < Rt^satj 


( 1 ) 


where Rosat = 0.13 is the saturation Rossby number, 
Rx = LxILboi, and Rx.sat = 10“^ '^ is the saturation Rx value. 
We use (3 - -2.7 (Wll). We assume that this relation can be 
used on t he PMS if the evolution o f Lboi and r* are treated 
correctly. ISanz-Forcada et aP (1201 ih derived a power law to 
convert Lx (5-lOOA) into EUV luminosity, Leuv (100-920A), of 
log Leuv = 4.8 -h 0.86 log Lx, where Lx and Leuv are in erg s“'. 

To calculate the evolution of the stellar radius, Lboi, the mo- 
ment of ine r tia, an d t*, we use the stellar evolution models of 
ISnada et^ (l2013h . However, their t* values are approximately 
a factor of two above those of W11 for 1 M© stars; we therefore 
normalize r* at all ages such that the MS value is consistent with 
Eqn.IU 


3. X-ray Observations 

To test our predictions for X-ray distributions at each age, 
we collect Lx values of single stars from ROSAT, XMM- 
Newton, and Chandra of open clusters with ages from 30 Myr 
to 620 Myr. The clusters are NGC 2547 (30 Myr ; iJeffries et al.l 
l2006h . a Persei (50 Myr; iProsser et al.l Il998h . NGC 2451 


(50 M yr; iHiinsch et al.ll2003h Blanco I (5 ( ) Myr ; |Pillitteri_et^ 
2003), Pleiades (100 Myr; [Micelaet^ 19991: j^auffer et ^ 


1994), NGC 2516 (110 Mvr: iPillitteri et al.ll2006l:Tleffries et al 


19971), NG C 6475 (300 My r: 


Prosser et al.lll995h . and Hyades 


(620 Mvr: IStern et al.iri995h. Tor NG C 6475, since no optical 
catalogue was available. iProsser et al.l (Il99-5h did not report up¬ 
per limits for the non-detected stars and therefore the percentiles 
for the distribution of Lx should be considered upper limits. 
Eor all MS clusters, except Blanco I where masses were given, 
we derive masses by conv e rting f rom (B - y)o using a relation 
derived from the I An et al.l (l2007h stellar evolution models. Eor 
the PMS cluster NGC 2547, we calculated masses using the 
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Fig. 2. Left (a): Predicted rotational evolution tracks for stars at the 10th (red), 50th {green), and 90th {blue) percentiles of the rotational distribu¬ 
tion. The solid and dotted lines show the envelope and core rotational evolution, respectively, and the horizontal solid lines show the observational 
constraints on the percentiles. The dashed b lack line shows the time dependent saturation threshold for M, B^ip, and Lx, calculated assuming a 
constant saturation Ro and the t* values of ISnada et al.l dTOLlh . Right (b): Predicted Lx along each of our rotation tracks and comparisons to 
observed LX values of single stars in several clusters, with upper limits shown by V symbols. The solid horizontal lines show the 10th, 50th, and 
90th percentiles of the observed distributions of Lx at each age, calculated by counting upper limits as detections. The two solar symbols at 4.5 Gyr 
show the range of Lx for the Sun over the course of the solar cycle. The scale on the right y-axis shows the associated L^uv- 


ISiess et ^ (l2000l) models. Since we use these X-ray observa¬ 
tions only to compare to our predictions from rotation, we do 
not attempt to homogenise the M* and Lx determinations for 
each cluster. Our quantitative determinations of the Lx tracks 
are based on the relation from W 11 where such homogenisation 
was done. 

4. Results 


(90th, blue), i.e., either as young PMS stars, as near-ZAMS stars, 
or as slightly evolved MS stars. The spread in Lx amounts to as 
much as 1.5 orders of magnitude for several 100 Myr. 

Fig. 3 gives the age when a star falls out of saturation, fsat, 
as a function of initial Q, derived from our rotation model. This 
“saturation time” can be approximated by 


= 2.9Q, 


1.14 
'0 ’ 


( 2 ) 


iJohnstone et al.l (1201 5l) combined rotation period measurements 
of four young clusters with ages of ~150 Myr and used a rota¬ 
tional evolution model to predict the evolution of the resulting 
distribution of Q on the MS. The sample contains 1556 stars in 
the 0.4-1.1 Mq mass range. In Fig. [T] we show predictions for 
the distributions of Lx based on these Q distributions at ages of 
150 Myr and 620 Myr comparing with observed values in the 
Pleiades and Hyades, respectively. There is good overall agree¬ 
ment, although intrinsic X-ray variability (typically factors of 2- 
3) introduces some additional scatter such as is visible for stars 
exceeding the saturation threshold. 

To predict the range of possible Lx evolution tracks, we cal¬ 
culate rotation models for solar mass stars at the 10th, 50th, and 
90th percentiles of the Q. distributions, shown in Fig. |2^. Our 
models fit well the observational constraints on the percentiles, 
except for a slight underestimation of the 10 th percentile track 
in the first 20 Myr. This might cause us to underestimate the age 
when stars on the 10 th percentile track come out of saturation 
by a few Myr. Fig. |2]r shows predicted tracks for Lx and Leuv 
together with observed Lx for stars in the 0.9 -1.1 Mq range 
for each cluster listed in Sect. |3] Due to the low number of ob¬ 
servations in NGC 2547 (30 Myr), we extend the mass range to 
0.8 - 1.2 Mq. The tracks correspond very well to the observed 
percentiles in the individual clusters given the somewhat limited 
observationa l samples. The solar Lx (6x10^^ - 5 x 1 0^^ erg s“', 
lAvresI (Il997h : iPeres et ^ (l2000l) : Lludge et aP (120031) :') has been 
included as well and fits our models excellently. 

Stars on our rotation tracks drop out of saturation at «i 6 Myr 
(10th percentile, red), «;20 Myr (50th, green), and a!300 Myr 


where fsat is in Myr and Qo is the rotation rate at 1 Myr in units 
of the solar rotation rate. Assuming that the saturation level. 
Lx,sat ~ 10 “^ *^Lboi,o 5 is constant in time, which is approximately 
true, we obtain log Lx,sat = 30.46. If we approximate Lx by a 
power law after fsat (see Fig. 2b), we obtain, for a given Qq, 


Lx - 


Tx.sat? if f ^ fsati 

at’’, if f > fsat- 


(3) 


We require that the power law also hts the Sun with 
erg s“’ at an age of f© 


Lx ,0 = 10^^-^ erg s ' at an age of f© = 4570 Myr. We thus find 


b ' = 0.35 logQo - 0.98, a - Lx.ef©". 


(4) 


For the 10th, 50th, and 90th percentiles in Qq, corresponding to 
Do ~ 1.8 Dq, 6 . 2 D 0 , and 45.6 Q© with fsat ~ 5.7 Myr, 23 Myr, 
and 226 Myr, respectively, we find 


Lx - 


2.0x 10^'r‘'2 
2.6 X 10^2f-1.42 
2.3 X 103®r2-5° 


^EUV 


7.4 X I0^‘r°®® 10th 

4.8 X 10^2f-i.22 5Qjjj 

1.2xl0“r2'5 90th 


where luminosities are in erg s '. The slope of the median Lx 
track, b = -1.42, is very close to values reported from lin¬ 
ear regression to the Sun in Time sample dGiidel et al.lll997[ 
iRibas et aDl2005h . These power-law hts, valid for f > fsat, thus 
describe the range of possible evolutionary tracks for Lx and 
Teuv- 
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Fig. 3. Saturation time as a function of initial rotation rate flo- To cal¬ 
culate each rotation track, we fit power laws to the core-envelope cou¬ 
pling timescale and disk-locking timescales for our 10th, 50th, and 90th 
percentile models of Tqe = 380^®'^"* and = 13.5f2g°'^, where the 
timescales are given in Myr and flo is in solar units. The dashed line 
shows our best fit, given by Eqn.|2l and the vertical lines show the satu¬ 
ration times of the 10th, 50th, and 90th percentile rotators. 

5. Discussion 

The large differences in the evolutionary tracks, and therefore 
Lx and Leuv values, make it necessary to reconsider critically 
the atmospheric erosion of planets by high-energy radiation. 
To first approximation, the thermal mass loss rate from a sim¬ 
ple hydrogen dominated planetary atmosph ere, Mni, can be es¬ 
timated using the en ergy limited approach dWatson et al.lllOSTI : 
iLammer et alJl2009t) . where Mpi is proportional to the incident 
stellar EUV flux for a given set of planetary parameters. We 
therefore assume that Mpi oc Feuv, where Feuv is the EUV 
flux at the planetary orbit. As an example, we consider the case 
of a 0.5 MEarth planet at 1 AU around a 1 Mq star with an 
initial hydrogen atmo sphere of 5 x 10“^ MEanh- For this case, 
iLammer et al.l d2014ll calculated Mpi from the atmosphere of 
3.5 X IO^^oth s“^ with Eeuv = 100 erg s“' cm“^, where is 
the mass of a hydrogen atom (see the fifth case in their Table 4). 
We therefore assume that Mpi = 5.9 x IO^Eeuv, where Eeuv is 
in erg s“^ cm“^ and M is in g s“'. 

We show in Eig.|4]the evolution of the planetary atmospheric 
mass between 10 Myr and 5 Gyr assuming that the central star 
follows the Leuv tracks shown in Eig. |2] Although in all three 
cases, the planet at 10 Myr has identical atmospheric masses, 
by 5 Gyr the atmospheric hydrogen contents are very different. 
Orbiting the slowly rotating star, the planet retains 45% of its 
initial atmosphere; orbiting the rapidly rotating star, the planet 
loses its entire atmosphere within 100 Myr; orbiting the 50th 
percentile rotator, the planet also loses its atmosphere, but this 
instead takes almost a Gyr. Although this is a simple calculation 
of a single example atmosphere, it is sufficient to show that the 
star’s initial rotation rate, and the subsequent rotational evolu¬ 
tion, is an important aspect that needs to be properly considered 
when studying the evolution of terrestrial planetary atmospheres. 
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Fig. 4. The evolution of the planetary atmospheric mass of a 
0.5 planet orbiting a 1 Mq star at 1 AU with an initial mass of 

5 X 10“^ Mganh- The tracks correspond to planets orbiting stars that are 
on the 10th (red), 50th (green), and 90th (blue) percentiles of the rota¬ 
tional distributions. The vertical lines show the stellar saturation times. 
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